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The Heisenberg antiferromagnet on a two-dimensional triangular lattice is a paradigmatic problem 
in frustrated magnetism. Even in the classical limit S — >■ oo, its properties are far from simple. 
The "120 degree" ground state favoured by the frustrated antiferromagnetic interactions contains 
a hidden chiral symmetry, and supports two distinct types of excitation. And famously, three 
distinct phases, including a colhnear one-third magnetisation plateau, are stabilised by thermal 
fluctuations in applied magnetic field. The questions of symmetry-breaking raised by this model are 
deep and subtle, and after more than thirty years of study, many of the details of its phase diagram 
remain surprisingly obscure. In this paper we use modern Monte Carlo simulation techniques to 
determine the finite-temperature phase diagram of the classical Heisenberg antiferromagnet on a 
triangular lattice in applied magnetic field. At low to intermediate values of magnetic field, we 
find evidence for a continuous phase transition from the paramagnet into the colhnear one-third 
magnetisation plateau, belonging to the three-state Potts universality class. We also find evidence for 
conventional Berezinskii-Kosterlitz-Thouless transitions from the one-third magnetisation plateau 
into the canted "Y-state", and into the 2:1 canted phase found at high fields. However, the phase 
transition from the paramagnet into the 2:1 canted phase, while continuous, does not appear to 
fall into any conventional universality class. We argue that this, like the chiral phase transition 
discussed in zero field, deserves further study as an interesting example of a finite-temperature 
phase transition with compound order-parameter symmetry. We comment on the relevance of these 
results for experiments on magnetic materials with a triangular lattice. 

PACS numbers: 67.80.kb, 75.10.-b, 75.10.Hk 



I. INTRODUCTION 

The problem of antiferromagnetism on a triangular lat- 
tice occupies a special place in the history of frustrated 
magnetism. The failure of the Ising antiferromagnet on 
a triangular lattice to order at any temperature has been 
widely celebrated since the pioneering work of Wannier^ 
and Husimi and Syoji^. Similarly, the hugely influen- 
tial idea of a zero-temperature quantum spin-liquid was 
first mooted by Anderson in the context of the spin- 
1/2 Heisenberg antiferromagnet on a triangular lattice'^. 
More recent calculations suggest that this model does or- 
der at T = 0, albeit with a much-reduced sublattice mag- 
netisation, in a coplanar three-sublattice "120-degree" 
state"*'^. However, even at a classical level, the finite- 
temperature physics of this ordered phase is far from 
simple. The 120-degree state possesses a (pseudo)-vector 
chirality defined by the handedness of the spin texture in 
each elementary triangle^. As a consequence, the model 
can support Z2 vortices as well as conventional spin- wave 
excitations, and its low-temperature phase has been ar- 
gued to be a "spin gel" in which both play an important 
role^^^'^. Despite a very determined effort in simulation, 
these ideas remain controversial^^"^^. 

Frustrated magnets also exhibit a fantastically rich 
range of phases in applied magnetic field, and once again. 



studies of the triangular lattice antiferromagnet have 
played a central role in forming opinion. The Ising an- 
tiferromagnet on a triangular lattice famously exhibits 
a one-third magnetisation plateau in applied magnetic 
field^^. A one-third magnetisation plateau is also found 
in the Heisenberg antiferromagnet on a triangular lat- 
tice, where it takes the form of a coUinear three-sublattice 
state stabilised by both thermal^'' and quantum^^ fluc- 
tuations. Fluctuations select two further phases as a 
function of magnetic field — a coplanar, three-sublattice 
"Y-state", which is a canted version of the 120-degree 
state, and a 2:1 canted phase, which is a coplanar, canted 
version of the one-third magnetisation plateau^^'^®. The 
same succession of phases also occurs in the XY anti- 
ferromagnet on a triangular lattice^^, and very similar 
magnetic phase diagrams occur in a wide range of other 
models. Consequently, the magnetisation process of the 
triangular-lattice antiferromagnet is often presented as 
the paradigm for the behaviour of a frustrated magnet 
under field. As such it serves as a useful starting point 
to discuss e.g., the classical Kagome^°'^^, and Shastry- 
Sutherland^^ antiferromagnets. 

Further motivation for studying triangular lattice an- 
tiferromagnets in field can be taken directly from exper- 
iment. A full magnetic phase diagram as a function of 
field and temperature has been measured for a range of 
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triangular lattice antiferromagnets^'^, which include the 
S=l/2 intermetallic GdPd2Al324, and the 5=5/2 insu- 
lating oxides RbFe(Mo04)2^^ and Rb4Mn(Mo04)32^. In 
all of these cases magnetic anisotropy^^ and/or weak in- 
terlayer coupling^^ must be taken into account. None 
the less, clear evidence is found in each case for the 
three phases found in a Heisenberg model — the Y- 
state, the one-third magnetisation plateau, and the 2:1 
canted phase. More general models, with competing 
or anisotropic exchange interactions, exhibit even more 
complex behaviour as a function of field and temperature. 
These remain a subject of intensive study for the exam- 
ples they provide of novel magnetic phases^^"'^^, and for 
their delicate and subtle finite-temperature phase transi- 
tions^^'^^. 

In the context of all this activity, it is perhaps sur- 
prising that there are only two published attempts to 
determine the phase diagram of classical Heisenberg an- 
tiferromagnet on a triangular lattice in applied magnetic 
field from Monte Carlo simulation^^'^^. And, while these 
authors agree as to the phases present, the nature of 
the phase transitions between them remains largely un- 
explored. In this paper we attempt to remedy this situa- 
tion by using modern Monte Carlo simulation techniques 
to study the phase transitions which occur in the anti- 
ferromagnetic Heisenberg model on a triangular lattice, 
as a function of temperature and magnetic field. The 
challenge — and interest — of this problem stems from 
the fact that the coplanar Y-state and 2:1 canted phases 
break both discrete symmetries of the lattice and spin- 
rotation symmetry in the plane perpendicular to the ap- 
plied field. Moreover, since the Mermin- Wagner theorem 
forbids the breaking of a continuous symmetry in two 
dimensions'"', spin-rotation symmetry is broken only at 
the level of a topological, Berezinskii-Kosterlitz-Thouless 
(BKT) phase transition^^''*^. This type of compound 
symmetry-breaking is notoriously difficult to disentangle 
in two dimensions. The example of coupled Ising [Z2] 
and XY [0(2)] fields, in particular, has a long history, 
dating back to work on spin glasses by Villain'''^. Ques- 
tions of ^2(8)0(2) symmetry breaking also arise in the XY 
antiferromagnet on a triangular lattice^^, and in models 
of coupled Josephson-j unction arrays^^. A central theme 
for each of these problems is when — if ever — Ising and 
XY symmetries are broken in a single, continuous, phase 
transition*^"^^ . 

In the case of the coplanar phases found in the 
triangular-lattice Heisenberg antiferromagnet in applied 
magnetic field, the relevant symmetry is Z3 ® 0(2), and 
we find that the compound nature of the order param- 
eters significantly modifies the phase transitions which 
separate coplanar phases from the high-temperature 
paramagnet. For low values of field, we find a dou- 
ble phase transition, with the system passing first from 
the paramagnet into the collinear one-third magnetisa- 
tion plateau, and then into the Y-state. These transi- 
tions belong to the three-state Potts and BKT univer- 
sality classes, respectively. Another BKT phase tran- 
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FIG. 1: (Color online) Magnetic phase diagram of the anti- 
ferromagnetic Heisenberg model on a triangular lattice, obtained 
from Monte Carlo simulation. Results have been extrapolated to 
the thermodynamic limit using finite-size scaling analysis, as de- 
scribed in the text. Continuous phase transitions are drawn with a 
dashed line, while Berezinskii-Kosterlitz-Thouless phase transitions 
are drawn with a dotted line. For fields ft < 3 a double transition 
is found upon cooling from the paramagnet, while for h > 3 only 
a single transition is found. Horizontal dashed lines represent cuts 
at fixed field, h = 0.4, h = 2.0 and h = 5.0, or fixed temperature 
T = 0.05, analysed below. The low-field region h < 0.2, left un- 
shaded, is beyond the scope of this work. The black open symbol on 
the T-axis marks the phase transition found at /i = 0, from Ref. 10. 

sition is found at roughly constant field, separating the 
one-third magnetisation plateau and the 2:1 canted state. 
Approaching saturation, a single continuous phase tran- 
sition is found from the paramagnet into the 2:1 canted 
state. This exhibits a non-universal jump in spin-stiffness 
and continuously varying exponents as a function of mag- 
netic field, and so does not belong to any conventional 
universality class. In reaching these conclusions, we pay 
careful attention to finite-size effects, which are found to 
be very large for low values of magnetic field. Our results 
are summarised in Fig. 1. 

The remainder of the paper is structured as follows : 
In Section II we briefly discuss the Monte Carlo method 
used and introduce the order parameters for the different 
phases, together with associated correlation functions, 
and finite-size ansazes. In Section HI we discuss the 
topology of the overall magnetic phase diagram, focusing 
on the importance of finite-size scaling. In Section IV 
we discuss the transitions between ordered phases at low 
temperature. In Section V we present representative 
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cases of double phase transitions upon cooling, charac- 
teristic of low and intermediate fields. In Section VI we 
discuss a representative case of the single phase transi- 
tion observed at high field and its properties. Lastly, we 
conclude with an overall summary of our results and a 
discussion of some of the remaining open questions in 
Section VII. 



II. MODEL, METHOD AND ORDER 
PARAMETERS 



The model we consider is defined by the Heisenberg 
Hamiltonian 



(1) 



where the sum runs over all nearest-neighbour 

bonds of a triangular lattice (assuming periodic bound- 
ary conditions) and the sum on i runs over N lattice 
sites. We consider antiferromagnetic exchange interac- 
tions J > 0, and spins are taken to be classical vectors 
of unit length. In the presence of magnetic field it is 
convenient to rewrite this Hamiltonian as 



A 



3J 



18J 



9J 
T 



h 

9J 



(2) 



where the sum on A runs over all triangular plaquettes 
in the lattice and 



N ^ ' 



(3) 



By inspection, for h < 9J the system takes on its mini- 
mum energy for 



h 

9J' 



(4) 



At zero temperature, this condition selects a manifold of 
three-sublattice states which evolves smoothly from the 
120-degree state at ft, = to saturation at h = 9J. In 
what follows we set J = 1, and measure h in units of J. 

The degeneracy of this manifold is lifted by thermal 
fluctuations, which select coUinear and coplanar states 
over the non-coplanar ones at low temperature, in a 
manifestation of the celebrated "order-from-disorder" ef- 
fect^^. The resulting phases are 

i) a coplanar Y-state with one spin pinned in the neg- 
ative direction and two canting "up" , i.e. a dis- 
torted version of the 120-degree state; 

ii) a collinear state at exactly h = 3, with two spins 
"up" and one spin "down" , i.e. a one-third magneti- 
sation plateau; 

iii) a coplanar 2:1 canted version of the plateau, which 
smoothly interpolates until the collinear saturated 
paramagnet is reached at h = 9. 



These phases are illustrated in Fig. 1. 



In the collinear one-third magnetisation plateau (ii), 
only the components of spin participate in symmetry 
breaking. In the presence of magnetic field, this phase 
breaks only a discrete C3 = Z3 symmetry of the lattice, 
and long-range order is permitted in 2D. However, the 
coplanar Y-state (i) and 2:1 canted phase (iii) also involve 
spin components in the S^-S^ plane. In this case, long- 
range order implies selecting a common plane for canting, 
breaking the 0(2) symmetry for rotation of spins about 
the direction of the magnetic field. Since the Mermin- 
Wagner theorem forbids the breaking of this continuous 
symmetry at any finite temperature in 2D, only a regime 
with "quasi-long-range" order, described by the algebraic 
decay of spin correlations, is permitted. Both coplanar 
phases also inherit the broken lattice symmetry, hence 
a compound (g) 0(2) symmetry is broken at T = 0. 
These phases with broken mixed symmetries can also be 
viewed as magnetic "supersolids"^^"^^. 



In order to study the finite-temperature properties of 
the model defined by Eq. (1) we perform large-scale 
parallel tempering^^ Monte Carlo simulations. Simu- 
lating this model is challenging because of the under- 
constrained nature of the T = ground state. At 
low temperature the system can become frozen in non- 
coplanar r = ground states which still obey the mag- 
netisation constraint, especially for small system sizes. 
In this case, the thermal-selection process is not fully 
realised and strong finite-size effects are visible. In or- 
der to overcome this problem, we couple the parallel- 
tempering Monte Carlo scheme to successive determin- 
istic over-relaxation sweeps, which comprise the reflec- 
tion of each spin around its respective local field. Since 
this is a reversible and micro-canonical update, the global 
Markov chain for parallel tempering and over-relaxation 
also obeys detailed balance on the whole. Simulations of 
from 48 to 128 replicas (temperatures) were performed in 
parallel for a variety oi L x L rhombohedral clusters with 
periodic boundary conditions. The linear size L was cho- 
sen to be commensurate with three-sublattice order in 
the range L S 60 — 210. Typical simulations involved 
2x10^ steps, half of which were discarded for thermalisa- 
tion. Each step consisted of one local-update sweep of the 
lattice followed by two over-relaxation sweeps, with repli- 
cas at different temperatures exchanged every 10 steps. 
Random initial configurations were employed. 



The three ordered phases found break the translational 
symmetry of the lattice. This can also be interpreted 
as the breaking of permutation symmetry between the 
three different sublattices in which the triangular lat- 
tice can be divided, labelled A, B and C. In order to 
study this process, we introduce a complex order parame- 
ter ^ = -\- i'ip2 , based on a two-dimensional irreducible 
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representation of the C3 ^ Z3 lattice rotation group 

ri = ^Y.^Sl + 2SI,~AS^c, (5) 



m' = \r,\' + \r2\\ 



A, (6) 
(7) 

where the sum over i runs over the N/3 elements of each 
sublattice A, B and C. Since parallel tempering effec- 
tively restores the lattice symmetries, we measure the 
magnitude of the order parameter 



(8) 



which is normalised to 12/\/Q in the case of a perfect 
"two-up, one-down" collinear configuration (only achiev- 
able at T = and h = 3). The ordering susceptibility 
and the (temperature-dependent) structure factor asso- 
ciated with this order parameter are defined as 



N 



(1^ 



T 



z\2\ 



T 



(9) 
(10) 



In order to characterise the phase transitions found we 
employ the standard finite-size scaling expressions 

O"- = L-'^/^O^'itL^/"), (11) 

X- = L^/''x"-(iLi/-), (12) 

S'-'' = L^-^'iS'^'itL^/"), (13) 

as a function of the reduced temperature t={Tc — T)/T. 
The critical exponents v, (3 and 7 are obtained from 
Eqs. (11) and (12), through the usual data collapse of 
the respective scaled quantities around the critical point. 
Since the scaled structure factor S^^ becomes indepen- 
dent of system size exactly at Tc, Eq. (13) allows the de- 
termination of the correlation-function exponent 77, pro- 
vided Tc is found beforehand. Error bars for the critical 
exponents are calculated by assuming a maximum devi- 
ation in the data points used to perform the data col- 
lapse, which results in a conservative (over) estimator of 
the precision in the obtained critical exponents. 

The presence of magnetic field reduces the symmetry 
of the model to S0(2) rotations in the S^-S^ spin plane, 
which supports, at most, quasi-long-range order. The 
spin stiffness ps, which acts as the (non-local) order pa- 
rameter for this algebraic order, can be calculated from 
the cost in the free energy of rotating the projection of 
each spin in the perpendicular plane = (S'f , Sf), see 
e.g. Ref. 34 and references therein. 



^/(ji:(e..,)S-xS 



(14) 



where Yij=Ti—Vj and ps has been normalised by 
\/3/2 to the unit area. Since parallel tempering ef- 
fectively restores the lattice symmetries, ps is aver- 
aged over three symmetric directions in the lattice 
e=(e„ ey)= {(1, 0), (1/2, ^3/2), (-1/2, ^3/2)}. 

The unbinding of vortex pairs at a BKT phase transi- 
tion suppresses the spin stiffness, leading to a jump in pg 
at the transition temperature Tbkt- This jump Aps can 
be expressed in terms of the correlation-length exponent 
r] as^^ 



Tf 



BKT 



2'Kri{Ti 



BKT) 



(15) 



In a conventional BKT transition, 77(Tbkt) = 1/4, and 
the jump in spin stiffness takes on the universal value 
Aps = 2Tbkt/t^- However, Tbkt is itself strongly renor- 
malised in finite-size simulations, and the correct ansatz 
for the finite-size scaling of Tbkt must take account of 
logarithmic corrections^^ 



Tr ^ T 



BKT 



1 

1 + 



1 



2 log L + log 6 



(16) 



Fits of this two-parameter scaling form to finite-size re- 
sults for the jump in spin stiffness give an estimate of the 
true Tbkt in the thermodynamic limit. 

The perpendicular component of the C3 order param- 
eter, can be defined by analogy with Eqs.(5)-(8). 
Since this quantity implies a (staggered) planar magneti- 
sation it must vanish at any finite temperature in the 
thermodynamic limit. However, in a BKT phase tran- 
sition, the scaling of the perpendicular structure factor 
, analogue to Eq. (10), yields the characteristic criti- 
cal exponent 77(Tbkt) = 1/4. 

Empirically, we find that the analysis of the correlation 
length ^ provides the most precise method to obtain the 
transition temperature for this model, being less sensi- 
tive to finite-size effects than, e.g.. Binder cumulants*^" . 
In order to calculate this quantity, the structure factor 
function associated with momentum q is first defined as 



Si exp(-iq.ri 



(17) 



Around the wave vectors corresponding to incipient 
three-sublattice order — Qk = {(47r/3,0), (27r/3, 7r/-\/3)} 
— and in a disordered phase, this quantity displays the 
characteristic Lorentzian form 5(q) oc ^^^^2 , arising 
from short-range correlations. For a sufficiently large 
system the correlation length ^ can be obtained from the 
ratio between S{q) at qk and the nearest allowed wave 
vector qK + Sq 



1 



5(qK) 



\Sq\ V 5(qK + Sq) 



1. 



(18) 



where Sq = (27r/T, 0). Eq. (18) is only directly related 
to the physical correlation length in the absence of long- 
range order, that is to say, of Bragg peaks in the structure 
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factor. The structure-factor function, and thereby the 
correlation length, can be divided into parallel (^^^), and 
perpendicular (^'^^) components, where the axis is 
defined by the direction of the magnetic field. 

The correlation length becomes infinite at any critical 
point, whether this is a conventional continuous phase 
transition, or the topological transition into a critical 
BKT phase. Hence the scaled quantities £,^^ /L or ^^^/i 
become independent of system size (or field) at this tem- 
perature, from where Tc (or he) can be found in an un- 
biased way. Error bars for are estimated using the 
difference between Tc as obtained by the intersection of 
the scaled ^ between the two largest clusters available 
and as obtained the intersection between the second- and 
third-largest clusters available. The data collapse of ^ 
for different system sizes can also be used to extract the 
correlation-length exponent v ai a, continuous phase tran- 
sition 

e = Lm"n- (19) 

III. TOPOLOGY OF THE PHASE DIAGRAM 

Our main results are summarised in Fig. 1. As men- 
tioned in the introduction, it is not our objective to ad- 
dress the properties of the peculiar Z2 iXiO(3) phase tran- 
sition at /i = 0, although it has been speculated that this 
survives the presence of a (very) small external field 
As we shall see, accurate simulations become increasingly 
difficult for very low fields and, owing to the very large 
correlation length in the S^-S^ spin plane, require sys- 
tem sizes larger than the ones presently available. We 
will therefore not discuss the phase diagram for h < 0.2, 
an area left unshaded in Fig. 1. 

For 0.2 ^ h < 3, two different phase transitions oc- 
cur as a function of temperature. At high temperature 
a continuous phase transition signals the breaking of the 
translational symmetry of the lattice along the spin 
direction. The resulting phase is the coUinear one-third 
magnetization plateau. We note in passing that the mag- 
netization of this plateau is not tied to one third in a clas- 
sical Heisenberg model at finite temperature, since the 
coUinear "up-up-down" state is dressed with thermally 
excited spin-wave excitations. (A perfectly coUinear "up- 
up-down" state is realised in a classical model only at 
T = and h = 3, where it is energetically degenerate 
with many other, non-coplanar, states). 

In the field range 0.2 < ft. < 3, the canted Y-state is 
found by lowering the temperature from the one-third 
magnetization plateau. The Y-state inherits the broken 
translational symmetry of the plateau and, since two of 
its spins are canted, also breaks spin-rotational symmetry 
at T = 0. At finite temperature, this results in a phase 
with algebraic order in the S^-S"^ plane. The best in- 
terpretation of our numerical results is that, in the ther- 
modynamic limit, the Y-state is never found to be in 
contact with the paramagnet, in agreement with other 
recent works'^ ^. 



For values of field above the plateau and below the 
saturation limit, 3 < /i < 9, a single continuous phase 
transition separates the paramagnet from the 2:1 canted 
state. This transition corresponds to the simultaneous 
onset of long-range order in the spin component and 
algebraic order in the S^-S"^ spin plane. 

The transition temperature associated with long-range 
order along the direction is obtained by the critical 
scaling of the corresponding correlation length The 
transition temperature between the plateau and the Y- 
state is obtained by the characteristic BKT finite-size 
scaling of the jump in the spin stiffness. Although the 
critical scaling of the correlation length should yield 
the same result, we find that this method is less accurate, 
even for intermediate fields. This can be explained by the 
very rapid growth of the /L ratio as field is lowered, 
until eventually this correlation length exceeds the linear 
sizes of the available clusters, rendering this analysis use- 
less. Hence, as we shall see, the determination of Tbkt 
using ^^'^ is less accurate than that using ps for the same 
set of system sizes, even at intermediate fields. 

The onset of algebraic order in the S^-S^ plane associ- 
ated with the 2:1 canted state is obtained through the on- 
set of critical scaling in both as a fixed-temperature 
scan from the plateau or as a fixed-field scan from the 
paramagnetic region. The transition temperatures ob- 
tained with the perpendicular (^^^ and parallel cor- 
relation lengths agree very well, although they relate to 
different symmetries. Moreover, a jump in the spin stiff- 
ness is observed very close to the transition temperature 
found with the analysis of the correlation lengths. 

These conclusions are valid in the thermodynamic limit 
L — 00. It is instructive to contrast them with the re- 
sults obtained for a finite-size cluster with L = 120, sum- 
marised in the pscudo phase-diagram Fig. 2. In order to 
determine the finite-size "transition" points we analyse 
the crossing of the correlation lengths between a single 
pair of system sizes, L = 108 and L = 120, and regis- 
ter the temperature where ps = 2T/tt for L — 120. At 
low field, the transitions between the paramagnet and 
the plateau, and between the plateau and Y-state, are 
now indistinguishable, within error bars. Hence the best 
interpretation is that a single phase transition separates 
the Y-state from the one-third magnetization plateau. 
In Fig. 2 the position of the plateau- Y-state transition, 
obtained with the critical scaling of f^^, is also shown. 
The difference between the transition temperature thus 
obtained and the one obtained using ps vanishes gradu- 
ally as field is increased. Together with the observation 
that the correlation length increases with decreas- 
ing field, eventually becoming larger than L — 120 for 
h < 0.8, we can understand that the determination of 
Tbkt from is very strongly affected by finite-size cor- 
rections. Therefore, we only plot the result obtained with 
Ps in the L — >■ c» phase diagram, cf. Fig. 1. 

When passing from the paramagnet to the 2:1 canted 
state at higher fields, a similar separation is observed 
between the transition temperature obtained using the 
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FIG. 2: (Color online) Finite-size pseudo phase-diagram of the 
antiferromagnetic Heisenberg model on a triangular lattice in mag- 
netic field, obtained from Monte Carlo simulation of an L X L rhom- 
bohedral cluster with L = 120. Continuous "phase transitions" are 
inferred from maxima in the correlation length along the direc- 
tion, and drawn with a dashed line. Berezinskii-Kosterlitz-Thouless 
"phase transitions" in the S^-S^ plane are obtained from either the 
position of the spin-stiffness universal jump for L=120, or from the 
scaling of the corresponding correlation length, and drawn with a 
dotted line. A double transition upon cooling is only observed at 
intermediate fields 0.6 J; ^ J; 3, contrary to the conclusions for 
L — > oo, cf. Fig. 1. A small difference between different measure- 
ments of the transitions in the S^-S^ plane is observed for both 
low and high fields. Purple dashed lines indicate cuts at fixed field 
and temperature analysed below. The low-field region, h < 0.2, is 
beyond the scope of this work. 



scaling of correlation lengths, and that obtained from the 
jump in spin stiffness. However, once again we observe 
that this difference vanishes when the correct finite-size 
scaling is performed. We return to this point below. 



set of algebraic order in the S^-S^ plane, as observed 
in the rise of the corresponding correlation length ^^'^ 
and spin stiffness, cf. Fig. 3(a) and (b). These simula- 
tion runs were performed without parallel tempering, but 
the over-relaxation procedure alone was enough to obtain 
good results. The transitions in this region are found to 
be rather easy to simulate, since correlation lengths are 
relatively small. 
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FIG. 3: (Color online) Double phase transition as a function of 
field, at a fixed temperature T = 0.05. The three different ordered 
phases, Y-state, m=l/3 plateau and 2:1 canted state, are found 
with increasing field. The magnetisation plateau is distinguished 
from the surrounding phases by the suppression of (a) the perpen- 
dicular correlation length and (b) the spin stiffness pg. Its 
structure is also revealed by the suppression of (c) magnetic sus- 
ceptibility X- Throughout this paper, lines connecting data points 
are guides to the eye, unless stated otherwise. 



IV. FIELD SWEEP AT LOW TEMPERATURE 

We start our analysis of the phase transitions with a 
scan in field at fixed temperature T = 0.05, cf. Fig. 3, 
which can also be found as a vertical dashed line in 
Fig. 1. All the ordered phases discussed are found as 
the value of field is increased : firstly we have the low- 
field Y-state, which is then followed by the one-third 
plateau aX h — 2.633(3), and finally the 2:1 canted state 
ai h = 3.010(3). Leaving the plateau by either lower- 
ing or increasing magnetic field corresponds to the on- 



The different nature of the magnetisation plateau, 
when compared to the surrounding canted phases, is 
clear in the suppression of magnetic susceptibility Xj cf. 
Fig. 3(c). This feature becomes more pronounced as the 
temperature at which the field scan is performed is low- 
ered. 

Both phase boundaries between the plateau and canted 
phase are approximately linear in temperature, and can 
be traced to different spin- wave excitations of the plateau 
state. The transition between the plateau and the 2:1 
canted state is only weakly dependent on temperature, 
indicating that both phases have roughly the same en- 



7 



tropy, arising from similar spin-wave excitations. The 
other spin-wave excitation inside the plateau corresponds 
to a canting of the "up" spins. Since this lowers the total 
magnetisation along S*^, this spin- wave is favoured ener- 
getically if field is decreased, making it more favourable 
to create a spin-wave than to "cant" all spins into the Y- 
state. This "protects" the plateau against the decrease 
of field, and it is this higher entropy of the coUinear 
phase^'' which makes the plateau- Y-state transition line 
slope downwards. 



V. LOW AND INTERMEDIATE FIELDS, h<3 




FIG. 4: (Color online) Double phase transition as a function of 
temperature from the paramagnet into the plateau and then into 
the Y-state, at h = 0.4. (a) The transition into the plateau breaks 
a Cs lattice rotation symmetry at a temperature Tc = 0.290(3), 
identified through the scaling of the correlation length [left 
inset] . The right inset shows the Binder cumulants associated with 
the C:j order parameter. The expected crossing at Tc is subject to 
very large finite-size effects, (b) The transition into the Y-state is 
heralded by the jump in spin stiffness ps, yielding Tbkt = 0.253(4). 
This value is obtained with a 1/logL scaling with system size, as 
shown in inset to (b). (c) The perpendicular correlation length 
is so large that no critical crossing is observed for the system 
sizes studied. All quantities are strongly renormalised by finite-size 
effects, even at low temperatures. 

For values of applied field in the range 0.2 < /i < 3 a 
double phase transition is found as temperature is low- 



ered from the paramagnetic region. In Fig. 4 we present 
results for h = 0.4, which is representative of the low- 
field region 0.2 ^ h < 1.2. The low-temperature region 
displays a finite value of the C3 order parameter O^^ 
[Eq. (8)], cf. Fig. 4(a), and spin stiffness ps [Eq. (14)], 
cf. Fig. 4(b). This is indicative of, respectively, long- 
range order in the direction parallel to field and quasi- 
long-range (algebraic) order in the perpendicular plane, 
as expected in the Y-state. However, the absolute value 
of these quantities at finite temperature is strongly de- 
pendent of the system size. This can be attributed to 
the proximity of the h = point, where both these or- 
der parameters vanish and the correlation length is very 
large (but probably finite)^". This explains the unusually 
strong finite-size corrections, which vanish gradually as 
field is increased. In fact, the most spectacular demon- 
stration of these problems lies in the absence of critical 
scaling of the correlation length, i.e. there is no com- 
mon crossing or collapse for different system sizes, cf. 
Fig. 4(c). Since there is no long-range order in the S^-S^ 
plane, Eq. (18) still provides an accurate estimation of 
the correlation length in that region (obviously, the 
same no longer holds for ^^^). The absence of a merger, 
or even a crossing, implies that the asymptotic regime, 
where the correlation length is infinite (i.e. ^^^/L w 1 
for finite clusters), has not been reached in the lattice 
sizes studied. This may be due to a slowdown in sim- 
ulation dynamics, arising from the pathological proper- 
ties of the h = point For higher values of magnetic 
field h « 0.8, a critical crossing is only observed in the 
largest pair of system sizes studied, L = 180 and L — 210. 
With increasing magnetic field this crossing is observed 
for gradually decreasing system sizes. 

These unusually strong finite-size effects make the ac- 
curate determination of the transition temperatures very 
hard. For h = 0.4, the critical scaling of the cor- 
relation length yields a reasonably well-converged value 
(i.e. the movement of the crossing point between succes- 
sive system sizes becomes smaller and is not significant 
for the largest ones employed) at Tc — 0.290(3), cf. left 
inset to Fig. 4(a). This phase transition corresponds to 
the breaking of translational symmetry. We reserve its 
characterisation to later in the text, for a value of field 
that allows a cleaner interpretation. The determination 
of Tc through the analysis of the Binder cumulants for the 

two-component order parameter, U = I — j^^jj^J^, 
converges noticeably slower to the thermodynamic limit, 
cf. right inset to Fig. 4(a). 

Analysing now the S^-S"^ spin-texture plane, the best 
fit to the evolution with system size of the position of 
the universal jump in spin stiffness is given by a log- 
arithmic form [Eq. (16)], with Tbkt = 0.253(4) and 
b = -1.7084(5), cf. inset to Fig. 4:(h). The resulting 
value for Tbkt is significantly lower than Tc, clearly im- 
plying an intermediate phase between the paramagnet 
and the Y-state that only breaks translational symmetry, 
i.e. the one-third magnetization plateau. A jump in spin 
stiffness at T sa 0.25 in the thermodynamic limit implies 
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a remarkable finite-size renormalisation of spin-stiffness, 
since ps is finite for tfie lattice sizes studied spin stiffness 
in a broad region above that temperature, cf. Fig. 4(b). 
This agrees with the strong variation of ps value with 
system size, as observed even deep inside the Y-state. It 
should be emphasised that the separation of these two 
transitions is observed only in the L — )■ oo limit. Taken 
at a fixed system size, both estimates for the transition 
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FIG. 5: (Color online) Double phase transition as a function of 
temperature from the paramagnet state into the coUinear m=l/3 
plateau and then into the coplanar Y-state, at /i = 2. (a) The 
transition into the plateau is heralded by the rise in the C3 order 
parameter. The transition temperature is = 0.364(1), obtained 
with the critical scaling of the correlation length in inset to (a), 
(b) The spin-stiffness pg jump signals the BKT transition into the 
Y-state at Tbkt = 0.138(3), a value obtained with a logarithmic 
scaling with system size, as shown in inset to (b). The paramagnet- 
plateau transition only displays a weak signature in (c) magnetic 
susceptibility, but shows a clear peak in (d) heat capacity. The 
position of the inner magnetic susceptibility peak in (c), scaling 
of the correlation length in inset to (c), and the inner heat 
capacity peak in (d) give a inaccurate estimative of Tbkt- 



temperature coincide at Tj. « 0.291(4) (within error bars) 
for the finite-size systems studied, cf. Fig. 2. It is tempt- 
ing, therefore, to infer that these two transitions appear 
to take place at the same temperature. However, since 
the correlation length is very large, any comparison at a 
fixed finite size is unreliable. 

A further manifestation of strong finite-size effects is 
seen in the presence of a small amount of non-coplanarity 
in the configurations found in equilibrium. This can be 
observed by looking at e.g. scalar chirality or quadrupo- 
lar spin moment quantities, which show a distinct sig- 
nature inside the y-phase [not shown]. These signals 
of non-coplanarity scale to zero with increasing system 
size, but very slowly. In spite of all these problems, the 
value of Tbkt extracted from the jump in spin stiffness 
Ps, obeys the logarithmic evolution with system size ex- 
pected for a BKT transition for fields as low as h — 0.2 
[see inset to Fig. 4(c)]. 

In order to obtain a clean characterisation of this dou- 
ble phase transition, we perform a similar analysis for 
a value of field h — 2, where these two transitions are 
now well separated, cf. Fig. 5. A sharp rise of the C3 
order parameter, cf. Fig. 5(a), is associated with the on- 
set of long-range order when entering the magnetisation 
plateau. This is observed to happen at T = 0.364(1), 
as found in the critical scaling of cf. the inset to 
Fig. 5(a). This phase transition can also be observed in 
thermodynamic signatures such as a very shallow sup- 
pression of the magnetic susceptibility [Fig. 5(c)], and a 
sharp peak in heat capacity [Fig. 5(d)]. 

A finite-size scaling analysis of this phase transition 
is performed using the critical exponents for the three- 
state Potts model in two dimensions v = 5/6 (correla- 
tion length), (3 — 1/9 (order parameter), 7 = 13/9 (or- 
der parameter susceptibility) and 77 = 4/15 (correlation 
function) [Ref. 61]. A perfect data collapse is obtained 
for the order parameter [Fig. 6(a)], order-parameter sus- 
ceptibility [Fig. 6(b)], and structure factor at T = Tc 
[Fig. 6(c)]. This unambiguously confirms that this phase 
transition belongs to the three-state Potts universality 
class for h = 2.0. 
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FIG. 6: (Color online) Data collapse for the continuous phase 
transition at ft = 2.0, T = 0.364, using the exact three-state Potts 
universality class critical exponents, (a) C3 order parameter, (b) 
C3 order parameter susceptibility and (c) associated structure fac- 
tor S"". 



At a lower temperature T « 0.15, the weak features in 
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heat capacity and magnetic susceptibility herald a tran- 
sition associated with quasi-long-range ordering in the 
5^-5^ plane. This can be interpreted as the formation 
of vortex pairs in the spin texture defined by the S^-S^ 
plane, arising from the two canted spins of the Y-state 
configuration. Therefore we find a BKT transition and 
respective rise in the spin stiffness, cf. Fig. 5(b). The 
transition temperature is found by tracking the position 
of the universal jump in the spin stiffness, Tbkt in inset 
to Fig. 5(b). Once again the best fit to the finite-size 
scaling is given by a logarithmic function of system size, 
yielding Tbkt = 0.f38(3) and b = 1.5051(4). 

In order to confirm the BKT character of this transi- 
tion, we study the critical scaling of 5^^, the structure 
factor associated with the component of the C3 order pa- 
rameter perpendicular to field. The crossing of 
for different system sizes at Tbkt (using the value found 
with the ps analysis) occurs for a value of the correla- 
tion exponent 77 = 0.26(2), cf. Fig. 7. This is in good 
agreement with 77(Tbkt) = f/4, the expected value in a 
standard BKT transition. 
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FIG. 7: (Color online) Evidence for Berezinskii-Kosterlitz- 
Thouless nature of the lower transition at /i = 2.0, T = 0.138. The 
collapse of the structure factor S^^ at the critical critical tempera- 
ture from Fig. 5(b) yields a critical exponent r] = 0.26(2), in good 
agreement with the BKT universality class. 

If we accept this finite-size scaling analysis at face 
value, we are lead to the conclusion that, for all val- 
ues of magnetic field 0.2 ^ ft. < 3, the I13 and 0{2) 
symmetries are broken at different temperatures. These 
two, distinct phase transitions are themselves perfectly 
conventional. The phase transition from the paramagnet 
to the collinear one-third magnetization plateau shows 
three-state Potts character, while the phase transition 
from the one-third magnetization plateau to the alge- 
braically correlated Y-state is BKT in nature. This is 
the interpretation given in the phase diagram Fig. 1. 

However this interpretation needs to be approached 
with some caution. Even for h = 2, some of the quantities 
calculated are strongly affected by finite-size corrections. 
This can be readily observed in the inset of Fig. 5(c), 
where it is shown that the value of Tbkt obtained from 
the critical scaling of the transverse correlation length 
depends on the size of the lattice studied. The absence of 
a good data collapse within the algebraically correlated 
Y-state is further evidence of strong finite-size effects. It 
would therefore be premature to rule out a single phase 



transition for /i — ?> 0. Given that the clusters used in 
the present simulations are not small (180x180 = 32,400 
sites), a fairly heroic act of simulation may be needed to 
finally resolve this question. 



VI. HIGH FIELD, h>3 

The behaviour of the Heisenberg antiferromagnet on 
a triangular lattice for magnetic field /i > 3 is dramati- 
cally different. The zero-temperature state is a 2:1 canted 
version of the one-third magnetization plateau, which 
breaks the translational symmetry of the lattice and spin- 
rotational symmetry in the — S*^ plane (at T = 0), 
i.c the same symmetries as the Y-state studied above. 
However, a single phase transition mediates between the 
high-temperature paramagnet and the 2:1 canted phase, 
in clear contrast with the case described in Section V. 

A selection of results for a representative field value 
ft = 5 is shown in Fig. 8. The rise of the C3 order param- 
eter, in Fig. 8(a), heralds the onset of long-range order 
along the direction. Strong finite-size artifacts are 
observed in the paramagnetic region close to the transi- 
tion, but the scaling of the correlation length yields a 
well-converged value of T^ = 0.212(1), cf. Fig. 8(b), with 
negligible finite-size corrections. The critical scaling of 
the perpendicular correlation length cf. Fig. 8(c), 
results in a value of Tj, = 0.212(1) (again with negligible 
finite-size corrections) which is in perfect agreement with 
the value obtained with We observe that this very 
good agreement, smaller than the statistical error bars 
AT = 0.001 - 0.004, is achieved for all values of fields 
ft > 3.3. More precisely, this agreement implies that the 
correlations along the direction and the S^-S^ plane 
become critical at the same temperature. This is good 
evidence for a single phase transition into the Y-state, 
without any intermediate phase. In the absence of any 
symmetry-breaking field in the S^-S^ plane, the rise of 
the respective correlation length corresponds to the on- 
set to algebraic order, cf. Fig. 8(c). This can also be 
observed in the rise of spin stiffness, cf. Fig. 8(d). 

Such a phase transition, which breaks a compound 
symmetry (8) 0(2), need not show the behaviour ex- 
pected of either a three-state Potts, or BKT transition. 
This phase transition retains a continuous character up 
to the largest cluster size studied L = 210, as can be veri- 
fied by the unimodal energy distribution at the calculated 
Tc, cf. inset to Fig. 8(a). We have also explicitly checked 
that the Binder cumulant for energy does not develop 
any characteristic signatures of a first-order transition as 
L — >■ cx) [not shown]. The heat capacity, cf. Fig. 8(e), 
only shows a broad peak at an higher temperature than 
the estimated Tc. Although the peak does not diverge 
with increasing system system size, the temperature of 
its maximum becomes lower. 

The strong finite-size corrections, observed at T > Tc 
in e.g. Fig. 8(a), make the finite-size scaling analy- 
sis much less precise than in Section V, but some con- 
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FIG. 8: (Color online) Continuous phase transition from the paramag- 
net into the 2:1 canted state as a function of temperature, for h — 5. (a) 
C3 order parameter, measuring the broken translational symmetry. In- 
set to (a) shows the single-peaked energy distribution at the transition 
temperature, (b) Scaling of the parallel correlation length 4^^, show- 
ing strong finite-size effects above the transition. Inset to (b) shows a 
detail of the scaling of which yields a crossing at Tc — 0.212(1). 
(c) The scaling of the perpendicular correlation length also shows 
strong finite size effects above Tc. This transition temperature is also 
found to be Tc — 0.212(1), as shown in the inset to (c). (d) The rise 
in spin stiffness ps also heralds the entry into the 2:1 canted state. In 
the inset to (d) a l/(27r77(Tc)) jump of ps/T is scaled as a function of 
1/L -\- resulting in — 0.217(4), in relatively good agreement 

with the scaling of the correlation lengths presented in insets to (b) and 
(c). (e) Heat capacity shows a broad peak at an higher temperature. 



elusions can still be reached. For h = 5, it is possi- 
ble to reliably extract the correlation-length exponent 
— 2.0(2), the order-parameter exponent /? — 0.50(5), 
and the correlation-function exponent 77 — 0.50(5) along 
the spin direction, cf. Fig. 9. This combination of crit- 
ical exponents does not appear to belong to any known 
universality class. 
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FIG. 9: (Color online) Scaling analysis of the continuous phase 
transition between the paramagnet and 2:1 canted state at /i = 
5,T = 0.212. (a) The collapse of simulation results for the 
correlation length yields the correlation length exponent u = 2.0(2). 
(b) The scaling of the respective component of the order parameter 
parallel gives u = 2.0(2) and /3 = 0.50(5). (c) The critical scaling 
of the parallel S^^ structure factor gives the correlation exponent 
r; = 0.50(5). 

The finite-size scaling analysis of the perpendicular 
component of the correlation length results in = 
2.0(3), cf. Fig. 10(a). This closely matches the v expo- 
nent obtained for the parallel component of the correla- 
tion length. This is in contrast with the BKT univer- 
sality class, where the correlation length diverges expo- 
nentially. The value obtained for the correlation function 
exponent r] = 0.27(2) at T = 0.212, our estimate for T^, 
cf. Fig. 10(b), is interestingly close to what is expected in 
a BKT transition. Nevertheless, this is, apparently, just 
a coincidence, as we shall see. 
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FIG. 10: (Color online) Scaling analysis of the continuous phase 
transition between the paramagnet and 2:1 canted state at /i = 
5,T = 0.212. (a) The collapse of the correlation length yields 
the correlation length exponent u = 2.0(2), exactly as in the S"^ 
direction, cf. Fig. 9(a). (b) The collapse of the iS^" structure factor 
gives the correlation exponent r] = 0.27(2). 



We anticipate that, regardless of the details of the 
phase transition, the low temperature 2:1 canted phase 
will contain bound pairs of vortices in the spin texture, 
which will unbind at T — T^, and precipitate a jump 
in the spin stiffness ps- However the fact that generi- 
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cally 77 7^ 1/4 means that the jump in the spin stiffness 
need not have the "universal" value Aps = ZT/tt [cf. 
Eq. (15)]. 

This interpretation is corroborated by the analysis of 
the data in Fig. 8(d) and its inset — the extrapolation 
to L — > 00 of the temperature at which the jump 
occurs yields the value Tc = 0.217(4). This value agrees 
within errors with the estimate Tc — 0.212(1) obtained 
with the scaling of the correlation lengths. The best fit 
to the finite-size data is now given by the power-law form 
a/L + b/L'^, not the l/log(L) scaling expected in a BKT 
transition. However if this transition does not belong to 
the BKT universality class, there is no good a priori rea- 
son to assume that Tc scales logarithmically with system 
size. 

The h = 5 results presented above are broadly repre- 
sentative of the phenomenology of the single transition 
from the paramagnet into the 2:1 canted state, in the 
range 3 < /i < 9. However, important details such as the 
critical exponents change as a function of field. 

h 



exponent a is negative for this phase transition. This ex- 
plains the absence of a singularity in specific heat at the 
transition temperature [cf Fig. 8(e)], and clearly distin- 
guishes the continuous phase transition here found from 
a three-state Potts transition. 

The variation with field of the spin correlations in the 
5^-5*^ plane at this transition are also noteworthy, cf. 
Fig. 12. The correlation-function exponent rj{Tc) is ob- 
served to also change continuously with field, albeit more 
slowly. We use this value of rj{Tc) to find the location of 
the non-universal jump in spin stiffness for the other val- 
ues of magnetic field at /i > 3.3 displayed in Fig. 1. The 
critical temperatures thus found are in good agreement 
with the ones obtained with the correlation length, tak- 
ing in account the uncertainties in determining rj. 
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FIG. 11: (Color online) Evolution of the critical exponents associ- 
ated with correlations of the 5^ components of spin at the continu- 
ous phase transition between the paramagnet and 2:1 canted state, 
as a function of magnetic field h. Both (a) the correlation-length 
exponent v and (b) the order-parameter exponent /3 increase with 
increasing magnetic field — lines are guides to the eye. However, 
(c) the ratio p/v is roughly constant at 0.24(3) — horizontal line. 

Both V and /3 exponents, associated with the com- 
ponent, increase monotonically with magnetic field, as 
can be observed in Fig. 11(a) and (b). However, the 
(3/v ratio remains roughly constant at f5/v — 0.24(3), cf. 
Fig. 11(c). The two-dimensional scaling law rj = 2(3/v 
therefore implies 77 ~ 0.5, which agrees with the result 
from the scaling of the structure factor at Tc for the 
studied fields, as exemplified in Fig. 9(c) for h = 5. The 
order-parameter susceptibility is strongly affected by 
finite-size effects (not shown). However, the critical scal- 
ing of yields a value of roughly 7 « 3 for /i = 5, 
which is consistent with the scaling relation 7 = 1/(2 — 77). 
The hyper-scaling law in two dimensions a = 2(1 — v), 
combined with Fig. 9(a), means that the heat-capacity 



FIG. 12: (Color online) Evolution of the correlation-function ex- 
ponent r)(Tc), related to the S^-S^ spin plane, of the single phase 
transition between the paramagnet and 2:1 canted state as a func- 
tion of field. The horizontal dashed line shows the standard BKT 
value r?(TBKT) = 0.25. 

Our results strongly suggest the existence of a point 
at T 0.31, h 3.2, where three apparently continuous 
transitions meet. Although simulations approaching this 
point become very difficult, we find no evidence for any 
of the transitions becoming first-order. 

Lastly, the zero-temperature phase transition from the 
(collinear) saturated paramagnet into the 2:1 state is 
observed at a field of ft. = 9. This corresponds to the 
opening of a gap, at the three-sublattice momenta {qk}, 
to spin-wave excitations inside the saturated phase. As 
the saturation field is approached the required computa- 
tional effort increases and the accuracy suffers. However, 
we also interestingly observe that the agreement between 
the Tc found by the scaling of both correlation lengths 
and Tc found by the spin-stiffness jump is better at both 
high and low values of the 3.3 < h < 9 range of applied 
field. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we have used modern Monte Carlo 
simulation techniques to explore the finite-temperature 
phase diagram of the classical Heisenberg antiferromag- 
net on a triangular lattice in applied magnetic field. The 
broad outline of this phase diagram has been known for 
some decades^^, and all of the phases predicted — a 
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classification 
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^ = 1/9 


T = 13/9 


1^ = 5/6 


,,=4/15 


n/a 


3-statc Potts 


p lat eau— Y- st at e 


2 


n/a 


n/a 


n/a 


n/a 


'J(^BKt) = 1/4 


BKT 


p.m. -2:1 canted 




/3 = 0.50(5) 


T»3 


■/ = 2.0(2) 


77=0.50(5) 


77(Tc) = 0.27(2) 


unlcnown 



TABLE I: Critical exponents for the different continuous phase 
transitions studied in this paper, as found from classical Monte 
Carlo simulation of the Heisenberg antiferromagnet on a triangular 
lattice, at different values of magnetic field h. 



coUinear one-third magnetization plateau, together with 
two algebraically-correlated coplanar phases, the Y-state 
and a 2:1 canted phase — have since been observed in 
experiment^''^"^^ . 

None the less, recent works not withstanding^^''^^, the 
nature of the finite-temperature phase transitions be- 
tween these phases remains surprisingly poorly under- 
stood. The interest of this problem lies in the fact that 
both the Y-state and 2:1 canted phase break two quali- 
tatively different symmetries — a discrete Z3 symmetry 
associated with their three-sublattice structure, and spin- 
rotation symmetry in S^-S^ plane. The order parameters 
for these phases therefore have a compound Z3 (g) 0(2) 
character, and it is interesting to ask how the two symme- 
tries are restored as the temperature is raised. Existing 
studies of this compound symmetry are very rare**®'^^. 
The answers we find, summarised in Table I, are remark- 
ably different for the two different phases. 

We consider first the case of low values of field h < 3. 
Here the system exhibits two phases — a long-range 
ordered one-third magnetization plateau, and an alge- 
braically correlated Y-state [cf. Fig. (1)]. Cooling from 
the paramagnet reveals two continuous phase transitions, 
the expected three-state Potts transition into the one- 
third plateau, and then a conventional BKT transition 
into the Y-state at lower temperature. Both of these 
transitions are well-characterised for /i = 2 [Fig. (5)]. 
However as /i — the correlation length in the S^-S'^ 
plane increases dramatically and, for h < 0.8, is com- 
parable with the linear dimension of the largest clusters 
simulated for all temperatures [Fig. 4(c)]. This effect, 
combined with the proximity between the three-state 
Potts and BKT phase transitions, makes the interpre- 
tation of simulation results extremely challenging. None 
the less, we are able to obtain a good finite-size scal- 
ing of results for spin-stiffness and correlation length 
down to /i«0.2, under the assumption that the two phase 
transitions remain distinct and well separated [Fig. 4(a) 
and (b)]. It is this, quantitative, analysis of the simula- 
tion results which leads to the phase boundaries shown 
in Fig. 1. 

At first sight, this result might seem to imply that the 
double phase transition — from paramagnet to plateau, 
and then from plateau to Y-state — survives all the way 
down to zero field {h = 0). We would however coun- 
cil caution : the physics of Heisenberg model in finite 
magnetic field may be very different from that in van- 
ishing field, where the order parameter and excitations 



have a qualitatively different character-'^" . We note that 
exchange anisotropy has recently been argued to act as 
a singular perturbation in the Heisenberg antiferromag- 
net on the triangular lattice^"^, and the same role may be 
played by applied magnetic field for the isotropic Heisen- 
berg model studied here. It could also be that proximity 
to the unconventional phase transition at /i = 0, ren- 
ders the finite size scaling used to extract phase bound- 
aries in Fig. 1 unreliable for /i — 0. Previous studies of 
the Heisenberg antiferromagnet on a triangular lattice for 
this range of fields have argued for both a single transi- 
tion from the paramagnet into the Y-state^^, and a dou- 
ble transition of the type described above, but with the 
temperature window between the two transitions closing 
as /i — )■ 0, see Ref. 21. Further simulations with larger 
cluster sizes, together with a more sophisticated analysis 
of results, will be needed to resolve this issue. 

Our results for higher fields, ft. > 3, point to a very 
different scenario. Here the system undergoes a single, 
continuous phase transition from the high-temperature 
paramagnet into the 2:1 canted phase. For this to hap- 
pen, the fluctuations in all three spin components must 
become critical at exactly the same temperature. This 
would not be unusual in a three dimensional frustrated 
magnet'^'*, but it has some very interesting consequences 
in the present, two-dimensional model. Considering first 
correlations of the components of spin, we find that 
the correlation length exponent v and order parame- 
ter exponent /3 increase with increasing magnetic field 
[Fig. 11(a) and (b)]. However, the correlation function 
exponent remains constant at w 0.5 [Fig. 11(c)], a value 
quite different from three-state Potts transition seen at 
low values of field [Fig. 6(c)]. Turning our attention to 
the — 5^ plane, the O^^ order parameter vanishes 
in the thermodynamic limit, and spin stiffness shows a 
jump at the transition temperature [Fig. 8(d)], as would 
be expected for a BKT transition. However in this case 
the correlation-length exponent, t]{Tc), varies with mag- 
netic field, and is generically different from the value 
?7(T'bkt) = 1/4 found at a BKT transition [Fig. 12]. This 
non- universal value of 77 implies a non-universal jump 
in the spin stiffness, and the temperature at which this 
jump in the spin stiffness occurs is found to scale as 
a polynomial in 1/L [inset to Fig. 8(d)]. These results 
point to a highly unusual line of continuous phase tran- 
sitions, interpolating from a point at which three critical 
lines meet {T « 0.31, h « 3.2), to the saturated state at 
(T = 0,ft = 9) [cf Fig. 1]. 

At an intuitive level, it is easy to see why a phase 
transition at which and 0(2) symmetries are broken 
simultaneously might be different from an isolated three- 
state Potts or BKT transition. The appeal to three-state 
Potts or BKT universality classes rests on the assumption 
of purely short-range interactions. This condition is un- 
likely to be met in the combined transition, where critical 
fiuctuations of one field can mediate a long range inter- 
action for the other. For example, vortices in the trans- 
verse components of spin carry a (topological) charge and 
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might be expected to couple to components of spins, 
invalidating the idea of short-range interactions between 
these Potts variables. And, conversely, these vortices can 
only exist inside a finite-size "box" set by the Potts de- 
gree of freedom. Long range dipolar interactions are well 
known to induce logarithmic corrections to scaling in con- 
ventional phase transitions^^, and it seems reasonable to 
suppose that long-range interactions modify the critical 
exponents in this compound phase transition. 

Precisely what happens where these lines of continuous 
phase transitions meet, at {T k, 0.31, h k, 3.2), is difficult 
to say, as it becomes increasing difficult to extract reli- 
able estimates of the critical exponents as this point is 
approached. Our best estimate of the exponent v as- 
sociated with correlations of the 5^ components of spin 
tends io V ^ \ as h ^ 3.2 from above [Fig. 11(a)]. This 
is roughly compatible with the three-state Potts value 
V = 5/6, seen for the transition from the paramagnet into 
the one-third magnetization plateau for /i = 2 [Fig. 6]. 
However, the order-parameter exponent tends to /? ~ 0.2 
[Fig. 11(b)], roughly double the three-state Potts value 
of j3 = 1/9. Similarly, considering spin correlations in 
the S^-S"^ plane, we find a correlation-function exponent 
ri{Tc) ~ 0.5 [Fig. 12], twice the value observed for the 
BKT transition from one-third magnetization plateau to 
Y-state for h = 2 [Fig. 7]. Clearly more work needs to 
be done to understand how the different phases come to- 
gether at this point. 

Conformal field theory (CFT) has proved to be a very 
powerful tool for understanding two-dimensional phase 
transitions®^. Within this approach, every continuous 
phase transition can be characterised in terms of a single 
parameter, the central charge c, and critical exponents 
are typically rational fractions with discrete values de- 
termined by c. In the present case, the three-state Potts 
transition from the paramagnet to the one-third magne- 
tization plateau for /i = 2 has central charge c = 4/5, 
while the BKT transition from the one-third magnetiza- 
tion plateau to the 2:1 canted phase for T — 0.05 has cen- 
tral charge c — 1. It has been argued that continuously- 
varying critical exponents arise most naturally in Gaus- 
sian (c = 1) CFT's with an additional, marginal, oper- 
ator^®. However it is hard to see how either a c = 1 
theory, or a direct product of a c = 1 theory with an- 
other CFT, can be reconciled with the variation of ex- 
ponents found in our simulations. We speculate that the 
transition from the paramagnet to the 2:1 canted state 
might therefore provide an example of varying critical 
exponents associated with a central charge c 7^ 1. This 
line of phase transitions also provides an example of the 
concept of "weak universality" , where the exponent i] 
is universal while /? and v are allowed to change. The 
idea of weak universality was first proposed in the con- 
text of the two-dimensional Ising model with four-spin 
interaction®'', which can be described by a c = 1 CFT®^. 
Our results suggest a further generalisation to these ideas 
to compound phase transitions which do not necessarily 
have c = 1 as the global charge. 



It is also instructive to compare these results with ex- 
isting work on Z2 (8) 0{2) phase transitions in two di- 
mensions. In principle, systems with Z2 0{2) symme- 
try breaking can also support a line of continuous phase 
transitions from the disordered state with continuously 
varying exponents'*®'*^. Recent numerical work suggest 
that this scenario is not realised in the most widely stud- 
ied model, the fully-frustrated XY model on a square 
lattice*^. However, in a recent twist to the story, Ising 
and BKT transitions have been observed to merge into 
a single phase continuous transition in a more general 
model^*'^^. 

The ultimate test of the results contained in this pa- 
per would be comparison with the magnetic phase dia- 
gram of a real triangular lattice antiferromagnet. Here 
the picture is obscured by terms not present in the 
isotropic Heisenberg model, notably magnetic anisotropy 
and coupling between triangular lattice layers^^. How- 
ever, published results for Heisenberg models with easy- 
axis anisotropy^'', easy-plane anisotropy®^, and interlayer 
coupling'^^ suggests that many of the most interesting fea- 
tures of the phase diagram Fig. 1 survive. Moreover, the 
rapid advances in experiments on cold atoms in optical 
lattices might make it possible to simulate a truly two- 
dimensional and isotropic Heisenberg antiferromagnet in 
the laboratory^". 

In conclusion, the behaviour of the Heisenberg model 
on a triangular lattice in applied magnetic field is much 
richer, and much less well understood, than usually sup- 
posed. In this paper we have used modern Monte Carlo 
simulation techniques to characterise the different phase 
transitions which occur as a function of temperature and 
magnetic field. The interest of this problems stems from 
the combined Z3 (g) 0{2) symmetry of low-temperature 
coplanar phases. For values of magnetic field h < 3.2, 
we find that the Z3 symmetry associated with three- 
sublattice structure and the 0{2) symmetry associated 
with the spin-rotations in the S^-S^ plane are broken 
at different temperatures. In contrast, for high values of 
magnetic field h > 3.2, we find that these symmetries are 
broken at the same temperature, in a line of continuous 
phase transitions with continuously varying exponents. 
Our results leave a number of important questions unan- 
swered, including the topology of the phase diagram for 
h ~^ 0, and the way in which Z3 and 0{2) symmetries 
combine for (T « 0.31, h « 3.2). Given the importance of 
finite size effects, it seems unlikely that these questions 
can be resolved by simulation alone, without further in- 
put from field theory. We therefore hope that this paper 
will help to re-open the discussion of this canonical prob- 
lem in frustrated magnetism 
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